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NOTATION 


P projection operator from X onto X 

X Banach space of periodic, continuous functions on [0,1] x [0,1] for the 
first case and similar definitions for the other cases 

X piecewise linear subspace of X 

X Banach space isomorphic to X 

Aif 5 A^A£; Ia^^I = l^ill^fl' 1^1 ~ l^ifl'^^ij* l^ifl " 

Aj:j S ijth rectangle of the domain 

(j) linear extensicj? of (po to X; <]i = ({'qP 

(j)Q a mapping creating an isomorphism between X and X 

o)e(6) modulus of continuity of the kernel function h(s,t) relative to s; 
oig(6) = sup|h(s + o,t) - h(s,t)|(s > 0, t i 1, la] i 6) 

V for all 

3 there exists 

II II norm in the appropriate space 

3 such that 

=> implies 

c’ (S) n times continuously differentiable on S 
p pressure 

u component velocity in x-direction 

V component velocity in y-direction 

A[ ] d^l ]/8 x 2 + a2[ ]/a2[ ]/ay2 

AA[ ] d^l ]/9 x 2 + 2 -^^4-+ 3‘*[ ]/9y‘* 

°x °y 

as boundary of domain S 

I 1. It ( 1 
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SUMMARY 


A general computer solution for the stationary Navier-Stokes equations 
is developed. It is explicitly proved that this solution converges to the 
true solution as the grid sige shrinks to zero. 


INTRODUCTION 


Existing algorithms for solving the Navier-Stokes equations liave one 
critical deficiency: if were the true solution and f the computer 

solution, tlve relationship between f and would be strictly heuristic 
[ ]. In this report, however, the relationship between \p* and | will be 
explicitly demonstrated, and it will be proved that | converges to ij)* as 
the grid size shrinks to zero. 


MAIN DEVELOPMENT 


The two-dimensional stationary Navier-Stokes equations are typically 
expressed by the set of equations 

uujj + vuy + qx ” + fi “ 0 (1) 

uVjj + vVy + qy - vAv + fa “ 0 (2) 

Ujj -f Vy “ 0 (3) 


with 


u(3s) ■ -bg , v(3s) “ bi 


( 4 ) 


Equations (1) and (2) express the conservation of momentum and equation (3) 
expresses conservation of mass, where u and v denote the velocity 


components in the x and y directions, respectively, and 8s denotes the 
boundary of the domain of consideration s. 

Solving the set of equations (l)-(4) is equivalent to solving the 
equation 


with 


vAAij) + ijiyAil^x + fly - ^2,5 “ 0 


i|)jj(8s) = bi , <|Jy(8s) » b2 


( 5 ) 


( 6 ) 


for the stream function <|j (see ref. 1), where -i|»y « u and = v. 

In reference 1, it is shown that solving equations (5) and (6) is equiv- 
alent to solving a sequence of linear partial differential equations: 




with 

t|/n,(3s) = 0 , Os) - 0 , m = 0,1,2, ... (8) 


where 


P(i|j) = vAAiJi + »|JyAli),j - t|jj^Aij<y + fi “ f2^ 


^ ~ 'P 

m m+1 ^m 


By inspection, we see that each element of this sequence is of the form 



wich 


$(98) - 0 , $n(9s) - 0 


( 10 ) 


Lemma Is 

A solution of equations (9) and (10) is a solution of equation (11) , and 
vice versa; 


Where 




f('l'm) - ~ f dS 

's 


( 11 ) 


and 


K,„(x’ ,y' ,x,y) « -AGy(x’ ,y ’ ,x,y)ij),,^Cx,y) + + 2Gxy('jim^^ “ 'Ji„i^^) 


'f'mjjy^^XX “ Gyy) 


and 


G denotes the biharnionic Green's function of the first type on S (see 
ref, 2 for the proof). 

Consequently, the initial problem “ equations (l)-(4) has been reduced 
to one of solving a sequence of Fredholm integral equations of the second 
kind, as expressed in equation (11). In particular, if ij)* denotes the true 
solution to equations (l)-(4), then 


lim tli^ “ - to 


Of course, one could stop with equations (9) and (10) and strictly solve 
them numerically, which has been done with great success for the nonstationary 
case in reference 3. However, if one were to computerize equations (9) and 
(10), one would end up in the following situation: 
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X, From reference 1, we would hove a sequence of solutions converg- 

ing to the true solution ij/*. Also, 


II !(,* - l|.„JI ' 2Z'"->” 


2 ""-' 


1 





2. If i|)mn were the actual computer solution for it would be 
practically impossible to determine the accuracy of relative to 

since equation (9) is an unbounded operator equation. 

However, because of the results of Lemma 1, the original problem has been 
reduced to one of solving a sequence of bounded linear operator equations; 
moreover, the equations are in a form in which it is possible to explicitly 
determine the accuracy between the computer solution and the true solution. 

The one real drawback with the integral equation representation is that the 
kernel is a function of the biharmonic Green’s function (see eq. (11)), an 
area primarily familiar to those in structural mechanics. 

Therefore, it remains to develop a computer solution with error estimates 
for equation (11). To accomplish this, we must first determine the character 
of the kernel of equation (11) . 

The kernel of equation (11) is of the form 

K_ “ + . . . 

hn y'* mjj 


In reference 4, it is demonstrated that G can be expressed in the form 


G 



log r + V 


with V a regular blharmonic function in S with continuous fourth deriva- 
tives throughout S. Let G = G + v, where 

G = r^ log r 
Sir 


By differentiating G, it follows directly that 


4 




yy 


Gxx 


(sin^ 0 - co8^ 0) 


and 


^xy 



cos 6 sin 


0 ) 


Therefore, the kernel K^CP.Q) of equation (11) is continuous except as 
P approaches Q and there it is of the form l/rpq» In reference 5, a 
numerical algorithm is developed for two-dimensional Fredholm integral equa- 
tions of the second kind with kernels of the form log(l/rpQ). For this ease 
it is proved (ref. 5) Chat the computer solution converges, as a function of 
grid st«;e to the true solution (as h > 0). It just so happens that the 
results of reference 5 also hold true if the kernel is of the form 1/rpq as 
P Q and continuous otlierwise. The argument goes through in exactly the 
same manner for l/r as it did for log(l/r) . In essence, by substituting 
X/r for log(l/r) in the development, and by slightly changing Che proofs 
for the new singularity, the results immediately follow. 

Let be the computer solution of — the solution of equa- 

tion (11) “• for partition A relative to reference 6. It therefore follows 
Chat 

ill ^ ill as A 0 

Also, relative to reference 1, as m where ijt’*' is the true 

solution of the problem (1-A) . Hence, 


$ni^ -K as 


A *>■ 0 , m « 


since 


II ij/„ * - i|)*|| S II - ijij^ll + II - '|)*ll ••►0 as m *>■ «> , A 0 
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Hej^ce, by applying the results proved In references 1, 2, 5. and 6, it 
has been possible to develop a computer solution aa a function of grid 
size A such that converges to \|j* as A approaches zero, providing 
of course that the conditions are satisfied as specified herein. 
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